Spatial coherence singularities of a vortex field in nonlinear media 
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We study the spatial coherence properties of optical vortices generated by partially incoherent 
light in self-focusing nonlinear media. We reveal the existence of phase singularities in the spatial 
coherence function of a vortex field that can characterize the stable propagation of vortices through 
nonlinear media. 



Vortices are known to occur in coherent systems having 
a zero intensity at the center and a well-defined phase 
associated with the circulation of momentum around the 
helix axis 0. The last decade has seen a resurgence 
of interest in the study of optical vortices 0, owing in 
part to readily available computer-generated holographic 
techniques for creating vortices in laser beams. 

If a vortex-carrying beam is partially incoherent, the 
phase front topology is not well denned, and statistics 
are required to quantify the phase. In the incoherent 
limit neither the helical phase nor the characteristic zero 
intensity at the vortex center can be observed. However, 
several recent studies have shed light on the question how 
phase singularities can be unveiled in incoherent light 
fields propagating in linear media |3,y). In particular, 
Palacios et al. [Jf used both experimental and numerical 
techniques to explore how a beam transmitted through 
a vortex phase mask changes as the transverse coherence 
length at the input of the mask varies. Assuming a quasi- 
monochromatic, statistically stationary light source and 
ignoring temporal coherence effects, they demonstrated 
that robust attributes of the vortex remain in the beam, 
most prominently in the form of a ring dislocation in the 
cross-correlation function. 

Propagating in nonlinear coherent systems, optical vor- 
tices become unstable when the nonlinear medium is self- 
focusing, and this effect has been observed experimen- 
tally in different nonlinear systems 0. However, sta- 
ble propagation of optical vortices in self-focusing non- 
linear media has been recently demonstrated, both the- 
oretically and experimentally, for the vortices created by 
self-trapping of partially incoherent light when the spa- 
tial incoherence of light exceeds a certain threshold j6] . 

The purpose of this Letter is twofold. First, we study 
numerically the spatial coherence function of a vortex 
beam propagating in a self-focusing nonlinear medium 
and reveal its importance for the study of singular beams 
in nonlinear media. Second, we provide a deeper physical 
insight into the effect of vortex stabilization by partially 
coherent light reported previously. 

We consider a phase singularity carried by a sin- 
gle vortex beam. To simulate numerically the prop- 
agation of partially coherent light through a nonlin- 
ear medium, we use the coherence density approach 



based on the fact that an incoherent light source can be 
thought of as a superposition of (infinitely) many coher- 
ent components Ej that are mutually incoherent, hav- 
ing slightly different propagation directions: E(r, t) — 
J2j Ej(r) exp(ikj_jr) exp(ijj(t)), where kj_j = k(aje x + 
PjBy) is the transverse wave vector of the j-th compo- 
nent, having direction cosines ctj and (3j, r = xe x + ye y , 
7j(t) is a random variable varying on the time scale of 
the coherence time of the light source, and k = 2ir/\ is 
the wavenumber. The vortex is introduced via a phase 
mask at the input face (z = 0) of the medium. To avoid 
complexities that may arise from incoherent light sources 
having abrupt boundaries, we assume the source has a 
Gaussian profile 
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Here 



9o is a parameter that controls the beam's coher- 
ence, A(r) — (r/wo) 2 exp(— r 2 /a 2 ) exp(i(p) is the vortex 
profile, and ip is the angular variable. 

Scaling the lengths in the transverse directions to 
xq = 1/im and the length in propagation direction to 
z Q = 2kx\, where we chose k = 27r/(230nm), the prop- 
agating field Ej(r,z) can be described by the nonlinear 
Schrodinger equation: 



■ dE 3 {v,z) 
dz 



V 2 x E j (r,z)+ri(r,z)E j (r ) z) = ) (2) 



where 7/(r, z) accounts for the nonlinear refractive index 
change in the material. We assume a photorefractive 
medium with a saturable nonlinearity having a response 
time much longer than the coherence time of the light 
source. In this case r\ depends on the time-integrated 
intensity. The intensity is given by I = J2j\Ej\ 2 . 
The nonlinear factor may thus be written in the form, 
T](r, z) — I(r, z)/(l + sl(r, z)), where s is a saturation pa- 
rameter. Whereas numerical solutions of Eq. @ may be 
readily computed using the coherence density approach, 
we shall later adopt the equivalent multi-mode theory p| 
to provide a physical basis for our findings. The numeri- 
cal simulations are performed using a split-step method. 
In the examples shown below we use N — 1681 com- 
ponents, wo = 1.8, a = 1.5 and s = 0.5. We chose 
0o = 0.64°. 
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Second-order coherence properties of the propagating 
field may be quantified by determining the the mutual 
coherence function r(ri,r 2 ;z) = (_E*(r 2 , z, t)E(ri, z, t)), 
where the brackets represent an average over the net field 
E(r, z, t) = Y^j=i -^i( r ' z ) ex P(*7j W)- Again, we assume 
that the random phase factors jj (t) vary on a much faster 
timescale than the response time of the medium. For 
the linear propagation, Palacios et al. Q| demonstrated 
that the phase singularities occur in the cross-correlation 
r(— r, r) of an incoherent vortex, where the origin of the 
coordinate system is chosen to coincide with the vortex 
center. 
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FIG. 1: Contour plots of the intensity (left column) and the 
modulus of the cross-correlation (right column) of an incoher- 
ent vortex with 9o = 0.64°. Contrary to the case of the linear 
propagation, there is a local intensity minimum in the beam's 
center. The cross-correlation, however, shows the same ring 
of phase singularities as predicted in the linear theory. The 
size is 35 x 35pim. 

In Fig. ^ we show the results for the nonlinear evo- 
lution of an incoherent vortex and the cross correlation. 
First, we notice that the beam intensity has a local min- 
imum in the center of the vortex, even after propagating 
many diffraction lengths. This is contrary to the case of 
the linear propagation where a beam with the same de- 
gree of coherence 9q has maximum intensity in the center 
of the vortex after only a few diffraction lengths. Also, 
if we had chosen to propagate an incoherent ring of light 
without topological charge instead of an incoherent vor- 
tex, we would also observe a maximum in the beam's 
center. Thus we can state that the coherence function 
of the vortex manifests itself in the intensity distribution 
of the light beam after propagating through a nonlinear 
medium. In fact, the intensity profile remains reminis- 
cent of a vortex, even if the intensity does not quite drop 
to zero in the center of the beam. 

Looking at the beam's cross-correlation, we clearly ob- 
serve a ring of phase singularities in the cross-correlation 
r(— r, r). Thus, as the first result of our numerical stud- 
ies we state that the phase singularities predicted for the 
incoherent vortices propagating in linear media also sur- 



vive the propagation through a nonlinear medium. This 
is not self-evident, considering that in the nonlinear case 
the single components that form an incoherent light beam 
do interact, contrary to the linear case. A physically in- 
tuitive explanation how this ring of phase singularities 
develops under linear propagation is given in Ref. Q. 
However, this issue becomes more complicated for the 
propagation in a nonlinear medium. 




FIG. 2: An incoherent vortex soliton calculated using the 
modal theory and three modes with the topological charges 
m — 0, 1 and 3. (a) shows the profiles of the three com- 
ponents, (b) shows the total intensity of the soliton and (c) 
shows its cross-correlation T(— r, r). 

Although the coherence density approach can be used 
to simulate the propagation of partially incoherent light 
with arbitrary precision, it is of little use when it comes to 
finding an explanation for the results obtained from the 
simulations. More physical insight into the problem can 
be obtained by using the modal theory of incoherent soli- 
tons 0- According to the modal theory, the incoherent 
solitons can be regarded as an incoherent superposition 
of guided modes of the waveguide induced by the total 
light intensity. Since the incoherent vortices that we are 
dealing with induce circularly symmetric waveguides, the 
guided modes we have to consider are also circularly sym- 
metric. To explain our numerical findings, we construct 
a partially incoherent vortex soliton using circularly sym- 
metric beams with topological charges m = 0, 1 and 2: 
^( r ) = Em=o E m(r) exp(irwp) exp(ij m (t)). This can be 
done using a standard relaxation technique [To| . A more 
precise modelling of incoherent solitons would require 
more modes. Here, we restrict ourselves to three modes 
only, assuming that for a partially incoherent vortex the 
m = 1 component should be dominant and that the next 
strongest components should be those with topological 
charge ml = m ± 1, i.e. m! = 0, 2. Indeed, we find that 
the main features of incoherent vortex solitons can be 
explained qualitatively using only these three modes. 

The relative intensity of the m — beam and the 
m = 2 vortex as compared to the m — 1 vortex con- 
trol the light coherence. However, in order to assure 
that the total topological charge of the light mtot = 
Im{(/ E*(r x V£')dr)}e 2 / J Idr is equal to one, we have 
to chose the m = and m = 2 components of equal in- 
tensity. In order to check whether this simple approach 
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FIG. 3: The intensity (left column) and the cross-correlation 
(right column) of the far field. The effects of the nonlinear- 
ity on the intensity distribution can be clearly seen, whereas 
the cross-correlation maintains more or less the structure one 
would expect in the case of linear propagation. 

yields the results that agree at least qualitatively with 
the full numerical model, we calculate an incoherent vor- 
tex soliton using the modal theory. The resulting shape 
of the vortex components, the total intensity, and cross- 
correlation r(— r,r) are shown in Fig. [21 Comparing 
Fig. ^ an d Fig- 121 we notice the presence of two simi- 
lar features: the local minimum of the intensity in the 
center of the beam, and the ring-like structure of the 
cross-correlation. Hence, these two phenomena can be 
explained by considering a simple modal representation 
of the incoherent vortex consisting of only three modes 
with the topological charges m = 0, 1 and 2. 

First, the local minimum in the center of the beam can 
be explained by the fact that the waveguide induced by 
the m = 1 and m = 2 components affects the m = 
mode in such a way, that it also develops a local in- 
tensity minimum in its center, a fact well known from 
the vortex-mode vector solitons. Second, the ring-like 



structure of the cross-correlation comes from the dif- 
ferent radial extent of the single components. As is 
known from the physics of vortex- mode vector solitons, 
the m = component has the smallest radial extent, 
whereas the m — 1 and m = 2 components have larger 
radii. Hence the cross-correlation given by T(— r, r) = 

^ft(-')^W) = El=o^(-^ m (r), is 
dominated for small r by the auto-correlated m = 
component, whereas the m = 1 component dominates 
for larger r. For even larger r, the m — 2 component 
can also come into play which can eventually result in a 
second ring of auto-correlation. 

Returning to the more precise, yet numerically more 
demanding coherence density approach, we show in Fig.|3| 
the situation in the far field. All parameters are identical 
to those used in Fig. ^ In the far field as well we ob- 
serve a ring-like structure of the cross-correlation func- 
tion r(— f,f), where f stands for the spatial coordinates 
in the far field. The intensity distribution in the far field 
can also show a local minimum in the center of the beam, 
contrary to what one would obtain if the vortex was prop- 
agating through a linear medium 4] , and also in contrast 
to the result we would obtain if we were propagating a 
light beam without topological charge. This emphasizes 
the importance of the interaction between the beam co- 
herence function and the nonlinearity. 

In conclusion, we have shown that the phase singulari- 
ties in the spatial coherence function found earlier in the 
linear optics can survive the propagation through nonlin- 
ear media when the singular beam creates an incoherent 
vortex soliton. Our results emphasize the importance of 
the spatial coherence function in the studies of the prop- 
agation of incoherent singular beams. Not only the phase 
structure, but also the intensity distribution strongly de- 
pends on the initial form of the coherence function of the 
light beam as it enters a nonlinear medium. 
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